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^ ' Abstract 

1—5 ! 

I ■ The evolution of polarized quark distribution functions is taken into account the 

T-H , gluon emission and absorption, quark pair production and annihilation processes 
and treated by a statistical method which provides quark distribution functions at 

■ small Q^. The evolution function of quark distributions is derived and solved by 

, numerical methods. Then the dependent and ag dependent structure functions 

^ I F2(x,Q^) are obtained and the example of it, the proton structure function, is 

\^ • calculated, whose increase and decrease tendency with Os at different x is consistent 

. with experiments. This is also an evidence that the statistical methods for the 

I nucleon is reasonable. 

p Following our statistical model proposed in [1, 2], we would like to extend our discus- 

^ ■ sion to the dependence of the quark distribution functions, but only one part of the 

^ . evolution equations in QCD. In the early discussion, we didn't take into account the de- 

• 1—1 ! tailed processes involving gluons, for example, gluon emission and absorption by quarks. 

^ I By a combination of effective field theory [1, 3, 4, 5] and statistical methods [2, 6, 7, 8, 9, 

d ' 10, 11, 12], the quark distribution at small is obtained as follows [2], 

, , 3M^Rx ff 1 dE 

q[x) = 



271 J M^exp(^) + 1^2' V J 

where, T, R, M are the nucleon's temperature, radius and mass, fi is the quark's potential. 

When is big, in the e-p deep inelastic scattering, the electron can probe many 
other detailed processes other than 'j*q q [13, 14, 15, 16, 17]. We assume that the 
dominant detailed process involves a photon and the total spin in the process is conserved. 
Therefore, for a quark with spin parallel to the total spin of the nucleon g |, we consider 
the following detailed processes: 

1. gluon emission: 7*g |— i> q ^ g, 

2. gluon absorption: q ^ g ^ q ^ 

3. pair production: j*g ^ q ^ q i, 

4. pair annihilation: g j g j— > 'j*g, 
while for g J,, they are: 

1. gluon emission: j*q q i g, 
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2. gluon absorption: q i g ^ q i 

3. pair production: j*g ^ q ^ q i, 

4. pair annihilation: q '] q ^*g. 

Thus the cross sections of the above mentioned processes should be added into the 
parton model cross section. So the modified quark distribution functions could be written 
as 

l2 



Ztt J X y y 

«s , f^dy X 

tts , Q"^ rdy ^ , ,y. 

ZTT /I^ J X X 

0!s , rdy y 



for g t, where is a cutoff value of the transverse momentum, and 

q i {x, Q')^qi {x) + ^log% f % i 

zTT IX J X y y 



as , rdy y 

-2^^°vio-^^^"^'^^^^^y 

, rdy y 



for q I, where g{x) is the gluon distribution function [2] (here we neglect the gluon 
evolution for approximation), 

_ leM'^Rx rf 1 dE 
3^"^'- 27r y M.exp(f)-i:B^- ^ ^ 

Pq^ql , Pqiqi are the splitting functions denoting the probability of emission of a gluon 
with a fraction of the momentum of the parent quark which emits it, and Pq^g, Pqig are 
the splitting functions denoting the probability of emission of a quark with a fraction 
of the momentum of the parent gluon. Thus the second to fifth items in cq (2), (3) 
denote the contribution of the gluon emission, pair production, gluon absorption and pair 
annihilation respectively. From the calculating of the cross section of the each process 
[18], we have 

and 

PqU{z)^Pqig{z)^^{z' + {l-zf). (6) 
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Substitute eq (1), (4), (5), (6) into eq (2), (3) we get: 

27rrfgT(x,Q2) 1 SM'^Rx 



as d log 



2g T (x) - -g{x) + J J T {z)dz, (7) 



where 

4 1 + 1 c/E _ 8 1 dE 

^ ^ ~ 3{l-z)z^J ^exp(^) + l^ 3{l-z)J M^exp(^) + 1^ 



and 



where 



3 ^2 y i^exp(f) 

as dlogQ^ 3 27r y X 

4 1 + ^2 M ^ 8 /-f 1 



fin- i— !— ^ / i _ 8 r 

^ 3(1 -z)z^ J M^f.:srr,(EzElL^ + ^ 3(1 -z) J i 



3 (1 - z)z^ J ^ exp(^) + 1 ^2 3(1 -z) J ^ exp(^) + 1 ^2 

z2 + (l_^)2.f 1 rf^ 



3 ^2 7 ^exp(f) -1 ^2' 
in which, the quark's potential could be written as [2] 

/^.u = ^K±N)> (9) 

where Vg is the valence quark number of flavor q and is the number difference between 
the quark with spin parallel and anti-parallel to the total spin of flavor q. 

Eq (7) and (8) cannot be integrated analytically except in the region x 1 or x ^ 0. 
At x=0.03, 0.06, 0.1, 0.2, 0.3 ,0.4 ,0.5 ,0.6, 0.7, Wc have calculated q(x) and g^jg^ 
numerically for each quark (u j, m l,u |, -u l,d ^,d l,d j, d |, s |, s l,s |, s j) inside a 
proton. We use the parameters that T = 63Mey, R = (ISOMeV)'^, M = 938MeV, 
A„ = 0.83, Ad = -0.43, A, = -0.10 [2, 19, 20, 21] for the proton. The results are shown 
in Table 1 and Table 2. 

From Table 1 and 2, we can get some conclusions: 

1. ~ '^diogQ^ quite the same at very small x for all quarks and at very large x for the 
quarks with small potentials. 

2. ^ 'd tog change from positive to negative at about x = 0.3 ~ 0.4 for dominant quarks 
(the quarks with big potentials) and at x around 0.5 for the quarks with small potentials. 

The structure function could be calculated by [18], [22] 



F2{x,Q^)^xY,e 



2 
1 

qqM 



^^^^ + 2^ V i '^^^ " 3^^^^ + J^'^^' 



(10) 
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where f{z) is defined by eq (7), (8). The running couphng constant as could be written 
as [23, 24, 25] 

47r 

^^ = —, Q^' (11) 

where is a renormahzation constant. Substitute eq (11) into eq (10), we get the 
dependent structure function: 



qqU 



(12) 

where as(/i^) is the couphng Constance at = fi^. 

We have plotted F2{x,Q'^) ~ x at = ^^(/i^), = 0.99as(/i^), = 0.9as(/i^), 
as = 0.5as(/i^) (see Figure 1, 2, 3, 4) and F2{x,Q^) ~ as at x = 0.03, x = 0.1, a; = 
0.2, X = 0.3, X = 0.4, X = 0.5, X = 0.6, x = 0.7 (see Figure 5) for the proton, which shows 
an agreement with the experiments [26]. 
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1: q{x) ~ X of each quark inside the proton (T = 63MeV, R — {180MeV) 
938MeV). 
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Table 2: g^g^ ~ x of each quark inside the proton (T = 63MeV, R = (180Mey)-S 
M = 938Mey). 
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Figure 2: F2{x, Q"^) ~ a; of the proton at ag — 0.99q;s(/u^). The unit of ag is ag 
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Figure 3: ^2(0:, Q^) ~ a: of the proton at ag — 0.9q;s(//^). The unit of as is q;s(/x^). 
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Figure 4: ^2(0:, Q^) ~ a: of the proton at ag — 0.5q;s(//^). The unit of as is ag 
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